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Abstract 

We study the production of bound hypernuclei 1 S-Be and 2 §-Mg via the (K~ , K + ) 
reaction on 12 C and 28 Si targets, respectively, within a covariant effective La- 
grangian model, employing E bound state spinors derived from the latest quark- 
meson coupling model as well as Dirac single particle wave functions. The if + S~ 
production vertex is described by excitation, propagation and decay of A and 
£ resonance states in the initial collision of a if ~ meson with a target proton 
in the incident channel. The parameters of the resonance vertices are fixed 
by describing the available data on total and differential cross sections for the 
p(K~ , if + )S~ reaction. We find that both the elementary and hypernuclear 
production cross sections are dominated by the contributions from the A(1520) 
intermediate resonant state. The 0° differential cross sections for the formation 
of simple s-state S~ particle-hole states peak at a beam momentum around 1.0 
GeV/c, with a value in excess of 1 fib. 

Keywords: Cascade hypernuclei, covariant model of (K~, K + ) reaction, 
quark-meson coupling model cascade spinors. 
PACS: 25.80.Nv, 24.85.+p, 13.75.Jz 



1. Introduction 

The study of the double strangeness (S) hypernuclei is of decisive importance 
for revealing the entire picture of strong interactions among octet baryons. The 
binding energies and widths of the S hypernuclear states are expected to de- 
termine the strength of the EN and SiV — > AA interactions, respectively. This 
basic information is key to testing the quark exchange aspect of the strong in- 
teraction because long range pion exchange plays essentially a very minor role 
in the S = —2 sector. Even though t-channel pion exchange between S and 
nucleon does operate, its strength is quite weak because the irEE coupling is 
smaller as compared to the itNN coupling |lj. This input is also vital for un- 
derstanding the multi-strange hadronic or quark matter. Since strange quarks 
are negatively charged they are preferred in charge neutral dense matter. Thus 
these studies are of crucial value for investigating the role of strangeness in the 
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equation of state at high density, as probed in the cores of neutron stars @, 0] 
and in high energy heavy ion collisions at relativistic heavy ion colliders (RHIC) 
at Brookhaven National laboratory 0, CERN [H] and FAIR facility at GSI @. 

The (K~ , K + ) reaction leads to the transfer of two units of both charge and 
strangeness to the target nucleus. Thus this reaction is one of the most promising 
ways of studying the S = — 2 systems such as S hypernuclei and a dibaryonic 
resonance (-ff),which is a near stable six-quark state with spin parity of + and 
isospin TjjS J3] • Several ways have been discussed to approach these systems in 
the past Many experimental groups have used the ( K~ K + ) reaction 

on nuclear targets to search for a H dibaryonic resonance 12, 3, EH EH, EH] • 

As far as S hypernuclei are concerned, there are some hints of their existence 



from emulsion events [17]. However, no S bound state was unambiguously 



observed in the few experiments performed involving the (K , K + ) reaction on 
a 12 C target 13, 1|| because of the limited statistics and detector resolution. 



However, in the near future experiments will be performed at the JPARC facility 
in Japan to observe the bound states of 5 hypernuclei via the (K~~ , K + ) reaction 
with the best energy resolution of a few MeV and with large statistics by using 
the newly constructed high-resolution spectrometers [l^|. The first series of 
experiments will be performed on a 12 C target. These measurements are of great 
significance because convincing evidence for the S single-particle bound states 
would yield vital information on 5 single particle potential and the effective SA 
interaction. Already, the analysis of the scarce emulsion [l9[ and spectrometer 
data [IH, EH have led to S-nuclear potentials with depths that differ by about 
10 MeV from each other. 

The (K~,K + ) reaction implants a S hyperon in the nucleus through the 
elementary process p(K~, A + )S~. The cross sections for the elementary reac- 
tion were measured in the 1960s and early 1970s using hydrogen bubble cham- 
bers I20j, l2ll, K22J, ba, uA |25j . The total cross-section data from these measure- 



ments are tabulated in Ref. [26j . In a recent study 27J , this reaction was inves- 
tigated within a single-channel effective Lagrangian model where contributions 
were included from the s-channel [see, Fig. 1(a)] and it-channel diagrams which 
have as intermediate states A and S hyperons together with eight of their three- 
and four-star resonances with masses up to 2.0 GeV [A(1405), A(1520), A(1670), 
A(1810), A(1890), 11(1385), £(1670) and £(1750), which are represented by A* 
and £* in Fig. la]. This reaction is a clean example of a process in which baryon 
exchange plays the dominant role and the i-channel meson exchanges are ab- 
sent, as no meson with S = +2 is known to exist. An important observation 
of that study is that the total cross section of the p(K~ , K + )Z~ reaction is 
dominated by the contributions from the A(1520) (with Lj2.j — -D03) resonance 
intermediate state through both s- and u-channel terms. The region for beam 
momentum (pk-) below 2.0 GeV/c was shown to get most contributions from 
the s-channel graphs - the it-channel terms are dominant only in the region p K - 
> 2.5 GeV. 

Almost all of the previous theoretical investigations of cascade hypernuclear 
production via (K~,K + ) reaction on target nuclei [H, [H, [2^, [30, 31 1 have 
used the framework of an impulse approximation where the hyperon production 
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(a) (b) 

Figure 1: (color online) Graphical representation of our model to describe p(K~ , K + )E,~ 
(Fig. la) and A(K~,K+) B -B reactions (Fig. 2b). In the latter case the shaded area depicts 
the optical model interactions in the incoming and outgoing channels. 



dynamics is separated from that of the relative motion in the entrance and 
outgoing channels. Thus, the hypernuclear production cross section is expressed 
as a product of the cross section of the elementary cascade production reaction 
and a term that accounts for the dynamics of the relative motion. None of these 
models has attempted to calculate the cross sections of the elementary reaction 
- they have been extracted from the sparse experimental data. Therefore, the 
results of these calculations carry over the ambiguities that are involved in the 
experimental values of the differential cross sections for the elementary reactions. 

In this paper, we investigate the production of cascade hypernuclei via the 
(K~ , K + ) reaction on nuclear targets within an effective Lagrangian model 32 . 
I33L I34J ] , which is similar to that used in Ref . [27[ to study the elementary pro- 
duction reaction, p(K~ , K + )E~ . We consider only the s-channel production 
diagrams (see Fig. lb) as we are interested in the region where pk- lies below 
2 GeV/c. The model retains the full field theoretic structure of the interaction 
vertices and treats baryons as Dirac particles. The initial state interaction of 
the incoming K~ with a bound target proton leads to excitation of intermediate 
A and £ resonant states, which propagate and subsequently decay into a S~ 
hyperon that gets captured into one of the nuclear orbits, while the other decay 
product, the K + goes out. In Ref. [27| . it was shown that six intermediate 
resonant states, A, A(1405), A(1520), A(1810), S, and £(1385), make the most 
significant contributions to the cross sections of the elementary process. There- 
fore, in our present study the amplitudes corresponding to these six resonant 
states have been considered. 



2. Formalism 

2.1. bound state spinors 

The S - bound state spinors have been calculated in the quark-meson cou- 
pling (QMC) model as well as in a phenomcnological model where they are 
obtained by solving the Dirac equation with scalar and vector fields having a 
Woods-Saxon (WS) radial form. In the latter case, with a set of radius and 
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Table 1: Parameters of the Dirac single particle potential (having a WS radial shape) for the 
S — bound and proton hole states. In each case radius (ro) and diffuseness (a) parameters 
were 0.983 fm, and 0.606 fm, respectively for both vector and scalar potentials. The binding 
energies (BEs) of the S _ states were taken from the predictions of the QMC model. The QMC 
BEs for the proton hole states are also shown together with the corresponding experimental 
values (given within the brackets). 



State 


BE 

(MeV) 


v v 

(MeV) 


V s 
(MeV) 


1 |-Be(ls 1/2 ) 
l-Mg(l Sl/2 ) 

2 l-Mg(lp 3 / 2 ) 

2 l-Mg(l Pl/2 ) 
12 C (lp 3 /2) 
28 Si(ld 5/2 ) 


5.681 

11.376 

5.490 

5.836 
14.329(15.957) 
10.071(11.585) 


118.082 
124.674 
167.124 
181.486 
382.598 
378.421 


-145.780 
-153.881 
-206.326 
-223.858 
-472.343 
-467.186 



diffuseness parameters, the depths of these fields are searched to reproduce the 
binding energy (BE) of a given state. Since the experimental values of the BEs 
for the S~ bound states are as yet unknown, we have adopted the corresponding 
QMC model predictions (as shown in table 1) in our search procedure for these 
states. Furthermore, the scalar and vector fields are assumed to have the same 
geometry. It should be noted that the depths of the potential fields in such a 
model are dependent on the adopted radius (ro) and diffuseness (a) parameters 
but there is no certain way of fixing them. Nevertheless, using same ro for all 
the states may make the search for the potential depths too restrictive. Some 
authors have used the root mean square radius (RMS) of a given state to fix 
the r parameter (see, e.g., Refs. [35| and [36j). However, such a procedure 
cannot be applied for the 2 bound states at this stage due to the lack of any 
experimental information about them. With these constraints, we show in Ta- 
ble 1 the resulting parameters associated with the scalar and vector fields of 
the phenomenological model for S~ bound and proton hole states for the two 
target nuclei. In this table the QMC predictions for the BE of the proton hole 
states are also shown. However, in the search procedure for these states the 
experimental values of the BEs (given within the brackets) have been used. 

The use of bound state spinors calculated within the QMC model provides 
an opportunity to investigate the role of the quark degrees of freedom in the 
cascade hypernuclear production, which has not been done in previous studies 
of this system. Since the cascade hypernuclear production involves large mo- 
mentum transfers ( 350 MeV/c - 600 MeV/c) to the target nucleus, it is a good 
case for examining such short distance effects. In the QMC model [37|], quarks 
within the non-overlapping nucleon bags (modeled using the MIT bag), interact 
self consistently with isoscalar-scalar (a) and isoscalar-vector (u>) mesons in the 
mean field approximation. The explicit treatment of the nucleon internal struc- 
ture re pre sents an important departure from quantum hadrodynamics (QHD) 
model 38]. The self-consistent response of the bound quarks to the mean a 
field leads to a new saturation mechanism for nuclear matter [37} • The QMC 
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model has been used to study the properties of finite nuclei [39( , the binding of 
u, rj, rj and D nuclei 4^, 41, 4^] and also the effect of the medium on and 
J/* production 

The most recent development of the quark-meson coupling model is the 
inclusion of the self-consistent effect of the mean scalar field on the familiar one- 
gluon exchange hyperfine interaction that in free space leads to the N — A and 



S— A mass splitting [44j. With this [45| the QMC model has been able to explain 
the properties of A hypernuclei for the s-states rather well, while the p- and d- 
states tend to underbind. It also leads to a very natural explanation of the small 
spin-orbit force in A-nucleus interaction. In this exploratory work, the bound 
E spinors are generated from this version of the QMC model and are used to 
calculate the cross sections of the 12 C(if-, K +) 1 S-Be and 28 Si(if , if+ ) 2 §-Mg 
reactions. 

To calculate the bound state spinors, we have used the latest version of 
the QMC model. In this version, while the quality of results for A and 3 is 
comparable that of the earlier QMC results 4l[, no bound states for the S 
states [45| are found. The latter is in agreement with the experimental obser- 
vations. This is facilitated by the extra repulsion associated with the increased 
one-gluon-exchange hyperfine interaction in medium. We refer to Ref. [45| for 
more details of this new version of the QMC. 

In order to calculate the properties of finite hypernuclei, we construct a sim- 
ple, relativistic shell model, with the nucleon core calculated in a combination 
of self-consistent scalar and vector mean fields. The Lagrangian density for a 
hypernuclear system in the QMC model is written as a sum of two terms, Cq\ ic 
= £qmc + C-qmci where 0, 

£qmc = i>N{r)[ii ■ d - M N (a) - (g u u(r) 

T N P 

+9 P \b{r) + -(1 + r?)A(r) ho]<Mr) 
+ i[(Vc(r)) 2 +m^(r) 2 ] 



+ ^[(V6(r)) 2 +m 2 6(r) 2 ] + i(VA(r)) 2 , ( I. ) 



and 



£qmc= ^ Y {r)[il-d-M Y {a)-{glu{r) 



-QMC — / j 

Y=A,S,H 

+g Y p llb{r) + eQ Y A(r) ) 7o ]^(r), (2) 

where ■0Ar('r), iPy(t), b(r) and u>(r) are, respectively, the nucleon, hyperon, the 
p meson and the to meson fields, while m a , m u and m p are the masses of the 
a, ui and p mesons. The A(r) is Coulomb field. g u and g p are the w-N and p-N 
coupling constants which are related to the corresponding (u,d)-quark-a>, 
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and (u, d) quark-p, g q pl coupling constants as g^ — 3g^ and g p = g q . I]f and Q Y 
are the third component of the hyperon isospin operator and its electric charge 
in units of the proton charge, e, respectively. 

The following set of equations of motion are obtained for the hypernuclear 
system from the Lagrangian density Eqs. U)-©: 

[*7- d- M N (a) - (g uU }(r) + g p ^-b(r) 

+|(l + 7f)A(r)) 70 ]^(r)=0, (3) 

[i 7 • 8 - M Y (a) - (g Y oj(r) + g p ljb{r) 
+eQ Y A(r)) 7o ]^ Y (r)=0, (4) 

(~V 2 r + ml)<r(r) = 

9 a C N {a)p s {r)+g Y C Y {a)p Y s {r), (5) 
(-V 2 . + ml)Lo(r) = g u p B {r) + g Y pl{r), (6) 

(-V 2 r + ml)b(r) = ^-p 3 (r) + gJlYp Y B (r), (7) 

(-V')A(r) = e Pp (r) + eQ Y pl(r), (8) 

where, p s (r) {p Y (r)), p B (r) (p^{r)), p 3 {r) and p p (r) are the scalar, baryon, 
third component of isovector, and proton densities at the position r in the 
hypernucleus • On the right hand side of Eq. ([5]) , a new and characteristic 
feature of QMC appears, arising from the internal structure of the nucleon 
and hyperon, namely, ^^^^(cr) = and g Y C Y (a) — — 9M q}' j ^ where 

g a = g a {a = 0) and g Y = g Y (u = 0). We use the nucleon and hyperon masses 
as parameterized in Ref . [45| . The scalar and vector fields as well as the spinors 
for hyperons and nucleons, can be obtained by solving these coupled equations 
self-consistently. 

In Figs. 2(a) and Figs 2(c), we compare the scalar and vector fields as cal- 
culated within the QMC model with those of the phenomenological model for 
lsi/2 S~ states of ^-Be and 2 |-Mg hypernuclei, respectively. It may be re- 
called that in the QMC model the scalar and vector fields are generated by the 
couplings of the a and u mesons to the quarks. Because of the different masses 
of these mesons and their couplings to the quark fields the scalar and vector 
fields acquire a different radial dependence. In contrast, the two fields have the 
same radial shapes in the phenomenological model. We notice that in general, 
the QMC scalar and vector fields are smaller in magnitude than those of the 
phenomenological model in the entire r-region. One interesting point to note is 
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Figure 2: (color online) [(a)] Vector and scalar potential fields for lsi/2 — state in Be. The 
QMC model and Dirac single particle results are shown by solid and dashed lines, respectively. 
[(b)] Moduli of the upper (|/|) and lower (\g\) components of the lsi/2 2 orbits in X =-B 
hypernucleus in momentum space. |/| and \g\ of the QMC model are shown by the solid and 
dashed lines, respectively while those of the phenomenological model by the dashed-dotted 
and dotted lines, respectively, [(c)] and [d] represent the same for 2 §--Mg hypernucleus. 



that for the heavier hypernucleus, both scalar and vector QMC fields have their 
maxima away from the point r = 0, in contrast to the phenomenological fields. 
In the mean field models of the finite nuclei the proton densities are somewhat 
pushed out as compared to those of the neutron, because of Coulomb repulsion. 
This causes the S~ potential to peak outside the center of the nucleus. This 
is a consequence of the self consistent procedure. In the case of a chargeless 
hyperon (e.g. A) such effects are not observed. 

In Figs. 2(b) and 2(d) the moduli of the upper and lower components of 
lsi/2 S~ momentum space QMC (solid and dashed line) and phenomenological 
(dashed-dotted and dotted) spinors are shown for the ^-Be and 2 |-Mg hyper- 
nuclei, respectively. It is seen that the spinors of the two models are similar 
to each other for momenta (fc) up to 2.0 fin -1 . Beyond this region, however, 
they start having differences. The position of minima in the phenomenological 
model spinors is shifted to higher values of k and their magnitudes are smaller 
than those of the QMC model. It should however, be remarked here that the 
structure of the minima reflects the size of the system. An improved search for 
the depths of the WS potentials in the phenomenological model as discussed, 
above might remove the differences seen between the spinors of the two mod- 
els. We further note that only for k values below 1.5 fm -1 , are the magnitudes 
of the lower components, \g(k)\, substantially smaller than those of the upper 
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components. In the region of k pertinent to the cascade hypernuclear produc- 
tion, \g(k)\ may not be negligible. Thus the relativistic effects resulting from 
the small component of bound states spinors could be large for the hypernuclear 
production reactions on nuclei (see also the discussions presented in Ref. (4^|). 

2. 2. Cross Sections for Hypernuclear Production 

In order to calculate the amplitudes (and hence the cross sections) of the 
hypernuclear production reaction (see Fig. lb), one requires the effective La- 
grangians at the meson-baryon-resonance vertices and the corresponding cou- 
pling constants, and also the propagators for various resonances. After having 
established these quantities the amplitudes of the graphs of the type shown in 
Fig. 1 can be written by following the well known Feynman diagrams and can 
be computed numerically. 

The effective Lagrangians for the resonance- kaon-baryon vertices for spin-i 
and spin- 1 resonances are taken as 

(1-v) 

£kbr 1/2 = -gKBB, 1/2 ^R 1/2 [xi^ l fK^ -^j—T'y fJ/ (d fl if K )]ip B , (9) 

£kbr 3/2 = 9KBRa/2 ^ d^xipB + h. c, (10) 

with M = (rriR ± mj), where the upper sign corresponds to an even-parity 
and the lower sign to an odd-parity resonance, and B represents either a nu- 
cleon or a S hyperon. The operator T is 75 (1) for an even- (odd-) parity 
resonance. The parameter \ controls the admixture of pseudoscalar and pseu- 
dovector components. The value of this parameter is taken to be 0.5 for the A* 
and E* states, but zero for A and £ states, implying pure pseudovector cou- 
plings for the corresponding vertices in agreement with Refs. E3] - It may be 
noted that the Lagrangian for spin-| as given by Eq. (10) corresponds to that 
of a pure Rarita-Swinger form which has been used in all previous calculations 
of the hypernuclear production reactions within a similar effective Lagrangian 



model |3al33|,l34| 



Similar to Ref. [13], we have used the following form factor at various vertices, 

A 4 + (s — m z y 

where m is the mass of the propagating particle and A is the cutoff parameter, 
which is taken to be 1.2 GeV everywhere which is the same as that used in 



Ref. [27 1 



The parameters of the resonance vertices were fixed in Ref. [27J by describ- 
ing the total cross section data on elementary reactions p(K~ : K + )'Er and 
p(K~, K )"EP, where the form of the spin-| interaction vertex was somewhat 
different form that given Eq. (10). In this paper, therefore, we recalculate the 
cross sections of the elementary reaction using the spin-| Lagrangian given by 
Eq. (10). Apart from the total cross sections, we also describe the differential 
cross sections of the p(K~, K + )Z~ reaction which was not done in Ref. [27j |. 
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Tabic 2: A and £ resonance intermediate states included in the calculations. 



Intermediate state 


Ll2J 


M 


Width 


9KRN 


9KR.S 


(R) 




(GeV) 


(GeV) 






A 




1.116 


0.0 


-16.750 


10.132 


£ 




1.189 


0.0 


5.580 


-13.50 


£(1385) 


Pis 


1.383 


0.036 


-8.22 


-8.220 


A(1405) 


Soi 


1.406 


0.050 


1.585 


-0.956 


A(1520) 




1.520 


0.016 


27.46 


-16.610 


A(1810) 


Pox 


1.810 


0.150 


2.800 


2.800 



The values of the vertex parameters were taken to be the same as those de- 
termined in Ref. [27j except for the vertices involving the £(1385) resonance, 
where the coupling constants (CCs) have been slightly increased in order to 
better describe the differential cross section data (see Table 2). 

The two interaction vertices of Fig. [1] are connected by a resonance prop- 
agator. For the spin-1/2 and spin-3/2 resonances the propagators are given 

by 



1(7^ + ^1/2) 

- _ iF /2)2 



Vr i/2 ~ ~2 TFZ 7^2' ( 12 ) 

and 



^-(m^-ir^/2)^ ' (16 > 
respectively. In Eq. (TT3]) we have defined 

P» v = 9^ - 5 W - iA—^V v + (pV - fT-f) ■ (14) 

In Eqs. (TT2"]) and (TT3]), Tr 1/2 and ^r 3/2 define the total widths of the cor- 
responding resonances. We have ignored any medium modification of the reso- 
nance widths while calculating the amplitudes of the hypernuclear production 
as information about them is scarce and uncertain. 

In the next section we describe the results of our calculations for the (K ~ , K + ) 
reaction on both proton and nuclear targets. 



3. Results and Discussions 

In Figs. 3a, we show comparisons of our calculations with the data for the 
total cross section of the p(K~ , K + )'E.~ reaction for K~ beam momenta (pk-) 
below 3.5 GeV/c, because the resonance picture is not suitable at momenta 
higher than this. It is clear that our model is able to describe well the beam 
momentum dependence of the total cross section data of the elementary reac- 
tions within statistical errors. The arrow in Fig. 3a shows the position of the 
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s-channel 

■ ■ ■ ■ u-channel 
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p k (GeV/c) 



C 80 
£ 60f 



( b ) p (K", K + ) S" 
1 .7 GeV/c 






60 r 


</) 




'qrl) 


40 t 


a 

TJ 


20 



TJ 


= 




-1 



-0.5 0.5 
COS(G) 



-0.5 0.5 
COS(G) 



(c) p (K\ K + ) E" 
,i 2.1 GeV/c 



Figure 3: (color online) (a) Comparison of the calculated total cross section for the 
p(K~ , if+)S — reaction as a function of incident K~ momentum with the corresponding 
experimental data. Also shown are the individual contributions of s- and u-channel diagrams 
to the total cross section. The arrow indicates the position of the threshold for this reaction, 
(b) and (c) Differential cross sections for the same reaction for K~ beam momenta of 1.7 
GeV/c and 2.1 GeV/c, respectively. 



threshold beam momentum for this reaction which is about 1.0 GeV. The mea- 
sured total cross section peaks in the region of 1.35-1.4 GeV/c which is well 
described by our model. We further note that the cross sections for p K - < 2.0 
GeV/c are dominated by the s-channel contributions. 

In Fig. 3b we compare our calculations with the differential cross section 
data of the p(K~, K + )E~ reaction for p K - values of 1.7 GeV/c and 2.1 GeV/c. 
These data were read from the corresponding figures given in Ref. [24]. Both 
calculated and experimental differential cross sections are normalized to the 
same total cross section. We see that our calculations describe the general 
trends of the angular distribution data well in the entire angular region for both 
the beam momenta. Nevertheless, a slight overestimate of the data is noted 
at the forward angles. There is a need to remeasure these differential cross 
sections at the JPARC facility to confirm and refine the old bubble chamber 



data of Ref. [24 1 



The beam momentum dependence of the 0° differential cross section (da/d£l)o° 
for the p(K~,K + )'Er reaction is an interesting quantity because it enters ex- 
plicitly into the expression for the cross sections of the (K~,K + ) reaction on 
nuclei (leading to the production of E hypernuclei) in the kind of model used 
in Ref. [19( . Hence, the beam energy dependence of the zero angle differential 
cross section of the hypernuclear production directly follows that of [(da/dSY) °]. 
In Fig. 4, we show the beam momentum dependence of this quantity (using the 
same normalization as those in Figs. 3b and 3c). We see that [(da/ dSi)o<>] peaks 
in the same region of px- as the total cross section shown in Fig. 3a. On the 
other hand, the situation regarding the momentum dependence of the available 
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1 1.5 2 2.5 3 
p K - (GeV/c) 



Figure 4: (color online) The zero degree differential cross section of the p(K , K + )"E reaction 
calculated as a function of beam momentum. 



experimental data on [(da / dSl)o°] is quite uncertain. The existing data reported 
in Refs. 21 , 2^, 24 1 differ considerably from each other. This may be due to nor- 
malization problems between different experiments or may be arising from the 
large errors in the Legendre coefficients. In fact only Ref. [24| shows the data 
explicitly for four values of p K - between 1.7 GeV/c to 2.64 GeV/c, together 
with the coefficients of the Legendre polynomial fits to the data. The other 
two references give only the coefficients of such a fit that have large correlated 
errors. The cross sections of Ref. [2l| could have maxima at both 1.4 GeV/c 
and 1.7 GeV/c within the statistical errors. The data of Ref. [23[ have a peak at 
1.74 GeV/c but the coefficients of the Legendre polynomial fits given there may 
have misprints - at one beam momentum they even give negative cross section. 
Therefore, the position of the peak in the experimental zero degree differential 
cross section of the p(K~, K + )'Er is uncertain. A proper measurement of this 
quantity at the JPARC facility would be very welcome in order to remove this 
anomaly. 

In calculations of the hypernuclear production reactions, we have employed 
pure single-particle-single-hole (EEp -1 ) wave functions to describe the nuclear 
structure part, ignoring any configuration mixing effects. The nuclear structure 
part is treated exactly in the same way as described in Ref. 33] . The amplitude 
involves the momentum space four component (spin space) spinors (ip) which 
represent the wave functions of the bound states of nucleon and hyperon. For 
the proton hole and S~ states, spinors generated within the QMC and the 
phcnomcnological models were used in the respective calculations. We have 
used a plane wave approximation to describe the relative motion of kaons in the 
incoming and outgoing channels. However, the distortion effects are partially 
accounted for by introducing reduction factors to the cross sections as described 
in Ref. 29]. Since our calculations are carried out all along in momentum space, 
they include all the nonlocalities in the production amplitudes that arise from 
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Figure 5: (color online) (a) Calculated differential cross section at 0° as a function of K~ beam 
momentum for the 12 C(K~ , K+) 12 ,- Be and 2s Si(K~ , K+) 2 %-Mg reactions. The spin-parity 
of the final hypernuclear states are indicated on each curve. The particle-hole configurations 
of the 1 |._Be(l — ) and the 2 =-Mg(2+) hypernuclear states are [(lf>3/2)p 1 1 (1*1/2)3— > 
and [(^d 5 ^2)p 1 > (l s i/2)s- ]2 + > respectively. For the 2 g_Mg(3 - ) state the results shown 
are the sum of the cross section corresponding to both [(1^5/2)^ 1 1 (1^3/2)3- an< ^ 
[(1^5/2)? 1 1 (lPl/2)s _ configurations. The solid and dashed lines represent the results 
obtained with the QMC and the phenomenological S — bound state spinors. Arrows show the 
threshold of the two reactions, (b) Contributions of individual resonance intermediate states 
as indicated near each line, to the zero angle differential cross section as function of K~ beam 
momentum for the 12 C(K~ , K + ) 12 ^- Be(l _ ) reaction. Their coherent sum is shown by the 
solid line. 



the resonance propagators. 

We have chosen the reactions 12 C(A'~, A+) 12 r-Be and 28 Si(K~ , A+) 2 1-Mg 
for the first application of our model. The reaction on the 12 C target is billed 
as the "day one" experiment at the JPARC facility. The thresholds for these 
reactions are about 0.761 GeV/c and 0.750 GeV/c, respectively and the momen- 
tum transfers involved at 0°, vary between 1.8 - 2.9 fm _1 . The initial states in 
both the cases are doubly closed systems. The QMC model predicts only one 
bound state for the S-Be hypernucleus with the S~ hyperon being in a lsi/2 
state with a binding energy as shown in Table 1. For the 2 §-Mg case however, 
three distinct bound S _ states with configurations lsi/2, lfa/2 and lpi/2 have 
been predicted. The binding energies of these states are shown in Table 1. It is 
evident that for this nucleus lj?3/2 and lpi/2 states are almost degenerate. This 
reflects the fact that the S-nucleus spin-orbit potential is weak. This is due to 
the fact that since the corresponding total potential depth is small, the gradient 
of this potential that contributes to the spin-orbit force is also small. 

In case of the 12 C target, the S~ hyperon in a lsi/2 state can populate 1~ 
and 2~ states of the hypernucleus corresponding to the particle-hole configura- 
tion [(1^3/2)^ \ (1 s i/2)h-]- The states populated for the 2 |-Mg hypernucleus 
are [2 + , 3 + ], [1~, 2~, 3~, 4~], and [2~, 3~] corresponding to the configura- 



12 



tions [(ld 5/ 2) p 1 ,(ls 1/2 ) H -], [(ld 5 /2) p \(lP3/2)H-], and [(ld 5 / 2 ) p \ (lp 1/2 )s-], 
respectively. In Fig. 5, we have shown results for populating the hypernu- 
clear state with maximum spin of natural parity for each configuration. The 
unnatural parity states are very weakly excited due to the vanishingly small 
spin-flip amplitudes for this reaction (see, e.g., Ref. for an extensive dis- 
cussion on this point). However, for the 2 |-Mg(3~) case, the results shown 
are the sum of the cross sections obtained with both [(ld 5 / 2 )~ 1 , (1^3/2)--] and 
[{^d 5 / 2 )~ 1 ,(lpi/ 2 ) s -] particle-hole configurations. The latter contributes sub- 
stantially (up to about 75% within our model) to the excitation of this state. 

In Fig. 5a, the 0° differential cross sections are shown as a function of the 
beam momentum that are obtained by using S~ bound state spinors calcu- 
lated within the QMC as well as the phenomenological model for the reactions 
12 C(K~ ,K+) 1 i-Be and 28 Si(if-, if + ) 2 l-Mg. The configurations of the final 
hypernuclear states are as described in the figure caption. In the calculations 
of our reaction amplitudes, the relative motions of K~ and K + mesons in the 
initial and final channels, respectively are described by plane waves. The dis- 
tortion effects, which primarily describe the absorption of the incoming K~ , 
are however, included by introducing factors that reduce the magnitudes of the 
cross sections. These factors are taken to be 2.8 and 5.0 for 12 C and 28 Si targets, 
respectively as suggested in Ref. 29]. This necessarily assumes that shapes of 
the angular distributions are not affected by the distortion effects. This aspect 
will be further investigated in a future study. 

We see that the QMC model cross sections are larger than those obtained 
by using the phenomenological model by about 10-15% in all the cases. This 
reflects the fact that in the region of momentum transfers relevant to these 
reactions both the upper and the lower components of the QMC spinors are 
higher in magnitude than the corresponding phenomenological ones. 

An important observation in Fig. 5a is that for both the hypernuclear pro- 
duction reactions, the cross sections peak at p K - around 1.0 GeV/c, which is 
about 0.25-0.26 GeV/c above the production thresholds of the two reactions. 
Interestingly, it is not too different from the case of the elementary S~ produc- 
tion reaction where the peaks of the total cross section as well as the zero degree 
differential cross section occur at about 0.35-0.40 GeV/c above the correspond- 
ing production threshold (see Figa. 3a and 4). Furthermore, the magnitudes of 
the cross sections near the peak position are in excess of 1 /ib. It is important in 
this context to note that the magnitude of our cross section for a 12 C target at 
a beam momentum of 1.6 GeV/c is similar to that obtained in Ref. (HI within 
an impulse approximation model. Moreover, our cross sections at 1.8 GeV/c 
also are very close those of Ref. 1191 for both the targets. However, we fail to 
corroborate the results of Ref. [19( where cross sections were shown to peak 
for p K - around 1.8 GeV/c. It is quite probable that the distortion effects are 
dependent on the beam momenta and may be relatively stronger at lower values 
of Pk- . Nevertheless, this is unlikely to lead to such a large shift in the peak 
position. In any case, this effect was not considered in Ref. [19] also. There may 
be a need to re-examine the beam momentum dependence of the zero degree 
differential cross section in order to understand this different. 
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In Fig 5b, we note that the contribution from the A(1520) intermediate state 
dominates the total cross sections over the entire regime of Pk- values. This is 
similar to that noted in the case of the elementary S~ production reaction. The 
A(1405), and £(1385) states make noticeable contributions only for p K - very 
close to the production threshold. Other resonances contribute very weakly. Of 
course, our results are quite dependent on the CCs of various vertices, which are 
somewhat uncertain. Nevertheless, the respective cross sections shown in this 
figure are robust. First of all these CCs provide a good description of the total 
as well a differential cross sections of the elementary S~ production reaction. 
Secondly there is very little scope for increasing further the individual contribu- 
tions of the A and E intermediate states, because the CCs of the corresponding 
vertices used by us are already larger than their upper limits suggested in the 
literature. Furthermore, the contributions of other resonances are too weak and 
even have the wrong p K - dependence. Therefore, the final results are unlikely 
to be affected too much by the known uncertainties in the corresponding CCs. 

4. Summary and Conclusions 

In summary, in this paper the cascade hypcrnuclear production reactions 
12 C(K~, AT+) 1 I-Be, and 28 Si(K~ , K+) 2 ^-Mg have been studied within an ef- 
fective Lagrangian model, using the proton hole and S~ bound state spinors 
derived from the latest quark-meson coupling model. This is for the first time 
that the quark degrees of freedom have been explicitly invoked in the descrip- 
tion of such reactions. We have considered the excitation of altogether six A 
and S hyperon resonance intermediate states in the initial collision of the K~ 
meson with a target proton. These states subsequently propagate and decay 
into a E~ hyperon and a K + meson. The hyperon gets captured in one of the 
nuclear orbits, while the meson goes out. We constrain the coupling constants 
at the resonance vertices by describing both the total and the differential cross 
sections of the elementary p(K + , K~)E~ reaction within a similar model. 

We have also performed calculations with the spinors obtained by solving 
the Dirac equation with vector and scalar potential fields having Woods-Saxon 
shapes (the phenomenological model). Their depths are fitted to the binding 
energies of the respective states (QMC model values for the S~ particle states 
and experimental values for the proton hole states) for a given set of geome- 
try parameters which are taken to be the same for the two fields. While for 
^-Be hypernucleus the shapes of the QMC fields are similar to those of the 
phenomenological model, the two differ considerably in the case of 2 §-Mg. For 
the cases studied in this paper, the hypernuclear production cross sections cal- 
culated with the QMC S~ spinors arc found to differ only slightly from those 
obtained within the phenomenological (the former being about 10-15% higher 
in magnitude than the later). The distortion effects are included by introducing 
reduction factors to the cross sections taken from the previous studies of this 
reaction. 

The zero degree differential cross sections for the S~ hypernuclear produc- 
tion reactions on the two targets considered here, have peaks around the beam 
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momentum of 1.0 GeV/c within both the QMC and the phenomenological mod- 
els. This peak momentum is above the corresponding production threshold by 
almost the same amount as the position of the maximum in the elementary to- 
tal as well as zero degree differential cross sections lies away from its respective 
threshold. The peak cross sections are in excess of 1 \xb. Furthermore, the to- 
tal hypernuclear production cross sections are dominated by the contributions 
from the A(1520) (.D03) resonance intermediate state which is similar to the 
case of the elementary S~ production reaction. Other resonances make notice- 
able contributions only at beam momenta close to the production threshold of 
the reaction. It is desirable to perform measurements for the differential cross 
sections of the elementary 'Br production reaction in a wide beam momentum 
range. 

This work was supported by the University of Adelaide and the Australian 
Research Council through grant FL0992247(AWT). 
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